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Abstract 

Analytical calculations have been used to describe the calorimeter event selection in small-angle electron- 
positron scattering at LEP1. The first-order QED correction to the Born cross-section has been derived 
with leading log and next-to-leading log approximations. The second- and third-order corrections are 
computed with leading accuracy. Our analytical results are illustrated by Tables and compared with the 
corresponding results obtained with the of Monte Carlo generator BHLUMI. 
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1 Introduction 



The test of the Standard Model (SM) based on the unified theory of elec- 
troweak interaction and quantum chromo dynamics, is one of the primary aims 
of experiments at LEP1 and LEP2. The small-angle Bhabha scattering (SABH) 
process is used to measure the luminosity of the corresponding electron-positron 
collider. Being the normalization factor, the cross-section for the SABH process 
affects all observable cross-sections and is a significant component in both the 
precision measurements of the SM parameters and investigation of new physics 
at LEP1 [1]. 

By now a purely experimental precision better than 0.1 percent for the lu- 
minosity determination has been achieved at LEP1 [2]. However, to obtain the 
total value a systematic theoretical error must be added. Ths latter is deter- 
mined by the accuracy of theoretical description of SABH cross-section, that 
takes into consideration the specidic features of event selection and the detector 
geometry at LEP1. Therefore, there is a need to decrease its precision to the 
level of experimental one or better it. It is no mere chance that lately much 
attention has been given to the theoretical investigation of the SABH process 
[3-14]. 

The theoretical calculation of the SABH cross-section at LEP1 involves two 
somewhat different tasks. The first one is to privide an adequate description of 
the experimental restrictions imposed on the event selection in terms of final- 
particle phase space. The second task consists in writing the matrix element 
squared within the required accuracy. There are two methods for theoretical 
investigation of SABH at LEP1: the method based on the Monte Carlo event 
generator [3-6] and the method resting on semianalytical calculations [7-14]. 

The advantage of the Monte Carlo method is that it enables simulation of de- 
tectors and event selection schemes of different types. The modern Monte Carlo 
generators employ some additional procedures (Yennie-Frautchi-Suura expo- 
nentiation [15], electron structure function method [16] and effective-coupling- 
constant approach [17]) to eliminate problems associated with the infrared di- 
vergence, and to take into account the leading corrections in the higher orders of 
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the perturbation theory. Although a part of the second-order next-to-leading 
correction is compensated due to the exponential form of the electron structure 
function (see, for example, [7]), the total second-order correction, including 
a 2 ln(Q 2 /m 2 ) terms (here Q 2 is the typical value of the transfer momentum 
squared), remains undeterminated. 

The advantage of the analytical method is that it allows the use of the exact 
matrix element squared, based on essential Feynman diagrams. The infrared 
problem in the context of this approach can be solved in the usual way by taking 
into account virtual, real (soft and hard) photon emission and pair production 
in every order of the perturbation theory. The shortcoming of the analytical 
approach is its low mobility with respect to the change in the experimental 
conditions used for the event selection. Nevertheless, the analytical calculations 
are very important because they privide checking many Monte Carlo calculations 
for different "ideal" detectors. 

At present, four types of "ideal" event selections are used for the comparison 
of various MC - generators. One of them is completely inclusive as to the final 
electron and positron. In this case, the event is the simultaneous detection of 
the scattered electron and positron with the invariant mass higher than a certain 
threshold value. The ring-shaped symmetric or asymmetric detectors are used 
for the detection, this method of event selection is called in Ref.[3] as BARE1. 

The other three methods of event selection, called in Ref.[3] as CALOl, 
CAL02 and SICAL2, are calorimetric. If the final particles (photons, elec- 
trons and positrons) move almost in parallel to each other, i.e., they form a 
cluster, then the total energy of all particles in the cluster is determined at the 
calorimeter event selection. If the final particles do not form a cluster, then 
CALOl and CAL02 detect only the scattered electron and positron, whereas 
SICAL2 does not discriminate photons and electrons (or positrons). CALOl 
and CAL02 differ in the shape and size of the cluster. The detectors of the 
same geometry as in the case of BARE1 are used particle detection. 

It is most natural to apply the analytical approach to the same methods 
of the event selection, which are used in the MC-generators. The analytical 
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formulae for the SABH cross-section at LEP1 were published for the BARE1 
method of event selection in the cases of symmetric [11-13] and asymmetric [14] 
detectors. Those formulae include the full first-order correction as well as the 
second-order correction with leading and next-to-leading accuracy. They also 
include the third- order correction in the leading approximation. It is just these 
contributions that must be calculated in order to achieve the required accuracy. 
Some brief analytical results for the CALOl and CAL02 event selections in the 
case of symmetric detectors have been given in Ref. [14,25]. 

In this paper we present full analytical calculations of the SABH cross-section 
at LEP1 for CALOl and CAL02 in the most general case of asymmetric de- 
tectors. The corresponding formulae include the same contributions as for the 
BARE1 method, except the second-order next-to-leading contribution. As to 
the SICAL2 method, we hope to investigate it in the following publications. 

The paper is organized as follows. In Section 2 we introduce the "observ- 
able" cross-section a b s with allowance for the experimental restrictions on the 
angles and energies of the detected particles. We also calculate the first-order 
correction for the BARE1 method in the case of asymmetric detectors. The 
formulae obtained are widely used below. In Section 3 the first-order correction 
is considered for the CALOl method, and in Section 4 - for the CAL02. In 
Section 5 we present the formulae for the leading contributions in the second 
and third orders using the electron structure function method. These formulae 
are universal because they are independent of the shape and size of the cluster. 
In Section 6 we present the Tables that illustrate the comparison of our analyt- 
ical results with the corresponding calculations made with the MC-generator 
BHLUMI [3]. 



2 First— order correction for BARE1 

For the investigation of radiative corrections to the Born cross-section of 
SABH at LEP1 it is convenient to introduce the dimensionless quantity 
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where Q\ = e 2 0\ (e is the beam energy, and 9\ is the minimal angle of the wide 
detector). The "observed" cross-section a b s is determined by the following 
formula in the BARE1 case 

f 7 7 ^ 7 2 j_ 7 2 ±r^cr^ c dcr(e + + e~ — > e + + e~ + X) 

CTo6 „ = / ds^erfydyera j^p^i — '- . (2) 

where X denotes undetected final particles, and cci (#2) and (f^ ((f^) are the 
energy fraction and the transverse component of the momentum of the electron 
(positron) in the final state. The functions 6? (i = 1,2) take into account the 
angular cuts, and the function takes into account the cutoff on the energies 
of the detected electron and positron 

e; = 0(0 3 -0-)0(0--0i), e c 2 = 0(0 4 - + )0(0 + - 2 ), e = 0{ Xl x 2 -x c ), 

0_ = ^, e + = ^. (3) 

X\E X 2 £ 

In the case of asymmetrical ring-shaped detectors we have 

03 > 04 > 02 > 01, Pi = % > 1 • 

Below we shall assume, for definiteness, that the final electron is registered by 
wide detector (03, 0i), and positron - by narrow one (04, 9 2 ). 

The Born contribution to E, for SABH process at LEP1, is determined as 
follows 

and in symmetrical case the limits of integration are (/o|,l). The Born cross- 
section is the same both for the inclusive and calorimeter event selections. For- 
mula (4) includes the contribution of the scattering diagram and its interference 
with the annihilation one. The contribution of the annihilation diagram is pro- 
portional to 0\ and may be omitted within the scope of 0.1% . When calculating 
the radiative corrections to the cross-section (4), we systematically omit the 
terms proportional to 9\. 

The first-order correction to the cross-section of electron-positron scattering 
includes the contributions connected with radiation of the virtual and real (soft 
and hard) photons and vacuum polarization due to the light fermions including 



hadrons. The contribution of vacuum polarization is universal, and therefore we 
shall not take it into account below. 



Ei — Y, v+ s + Eff + E 



H 



(5) 



The correction due to the virtual and soft real photon emission (with the energy 
less than As, A <C 1) can be written as follows [12] (in this case X\ = x 2 = 

i, <?T + <?i = o) 

3. ^ , , e 2 d\z 



2 

o P r dz 



^v+s = 2- / -^[2(L-l)ZnA + ^L-2], L = In 



(6) 



where z = q^/Qi, m is the electron mass. 

The second term on the right side of Eq.(5) is due to the emission of the 
hard( with energy greater than As) photon by the positron registered by the 
narrow detector. In this case we have 



X = 7(/c ± , 1 - x 2 ), xi = 1 , k ± + qf + qf = 0, x c < x 2 < 1 - A 



(7) 



where k is the transverse component of the three-momentum of the hard pho- 
ton. The corresponding contribution can be obtained by the integration of the 
bremsstrahlung differential cross-section in the electron-positron collision over 
the region 



?J-2 



1 < z < p\, x 2 pl < z\ = ( y- 2 - < x 2 p\, — 1 < cosip < 1 

Qi 



(8) 



where if is the angle between two-dimensional vectors q± and q 2 . The result is 
defined by the formula 



a 



2tt j x z~ 



p3 dz ] - A i 



1 + x* 
1 — X 



dx 



(L - 1)(A 42 + A { £) + If (a, z; p 4j p 2 ) 



(9) 



where 



= ^-4^(A 42 + Ag) + a 42 L, + AgZ 2 + - e { 2 x) )L 3 + (# 4 - # 2 )Z 4 



1 + £ z 



Li = In 
Z 3 = In 



(x/9 4 — Z){xp\ — Z) 
[z — X 2 p 2 ){xp\ — Z) 



, L 2 = In 



(^-X 2 p 2 )(^ P4~^) 



{x 2 p\ — z){xp\ — z) 



La = In 



x 2 (xp\ — z){xp\ — z 
(z - p\){xpl - z) 



(pi - z){xp\ - z) 



(10) 



and the following designations are used for 9 - functions: 

£s$ = oP-o$\ A A2 = e 4 -e 2 , 
e\ x) = e(x 2 P 2 i -z), e i = e( P 2 -z), m^ = i-e\ x \ e i = \-e i . (11) 

The term E H , related to the emission of the hard photon by electron registered 
by means of the wide detector, can be obtained from (9) by change of the 
integration limits over z: pf) instead of (p 2 , 1) and substitution 

A42-1, K(x,z-p 4 ,p 2 )^K(x,z-p 3 ,l) (12) 

in the integrand. It is easy to see that the contribution, including the analogue 
of L4 , vanishes, and it is caused by the specific character of the accompanying 
9 - functions. 

The individual parts on the right side of Eq.(5) depend on the auxiliary 
parameter A, but this dependence disappeares in the sum, so that the first- 
order correction can be written as follows 

a ( P r dz \ r 1 1 \ 1 4- r 2 

Si = -{/^ / [(--^(l-x)+(L-l)P 1 (x))(A 42 +Ag)+ T ^if (x, ^; P4, P2) 

1 %c 

+ 1%) f(-kl-*) + (£-l)^ 

J „ z z \ L 2 1 — x J > 

(13) 

where P\{x) determines the iterative form of the nonsinglet electron structure 
function (see, for examle, Ref.[7]) 

P^x) = \^~0(1 - X -A) + (2/nA + hs(l - x), A -> . 

J. 00 Zi 

In order to make the cancellation of the A - dependence on the right side of 
Eq.(lO) more transparent, one may use the following relations 

1 x c -1 1 2 1 y/z/ps-t _|_ 2 

J P 1 (x)dx = - J — — -<fe , J P 1 (x)6 ( i ;) dx = 6 ( z c) J —^-dx, 

— X 

Jy C [J O-C ^ C 

1 y/z/pi-. . 2 Vz/p2-t 1 2 

/ p 1 (x)A ( $dx = e A e { t ) J —^-dx-d^ J —^-dx, 

rf, >• I 00 rf, l 00 

JjQ OjQ 

^42 — ^42 — ^42 • 



dx 



As we already noted, the right side of formula (13) defines the first-order cor- 
rection to the Born cross-section of SABH at LEP1 with switched off vacuum po- 
larization. The last one can be included by inserting the quantity [1 — U(zQl)]~ 2 
in the integrand (as to II see Ref. [3] and the bibliography cited there). 

The upper line on the right side of Eq.(13) corresponds to the emission of 
the real and virtual photons by a positron and lower one - by an electron. 
Besides, the terms that are accompanied by the ^-depending (rr-independing) 
6 - functions describe the emission in the initial (final) state. This information 
is very important for the investigation of calorimeter event selection. 

But before to move on this investigation , we note one more circumstance. 
When deriving the Eq.(13), we systematically ignore the terms of the order of 
Q 2 « Q 2 /S (S = 4e) in comparison with a unity. However, as it is known [18], the 
terms of this type have a double-logarithmic asymptotics. Their contribution 
to Si equals parametrically to [12,19] 

aQ 2 2 Q 2 

and it amounts to about 10~ 4 in the LEP1 conditions. Therefore, we expect 
the systematic deviation of our calculations in the first order of the perturba- 
tion theory from the corresponding results of MC-generator BHLUMI, which 
completely includes such contributions on the 0.01% level (see Section 6). 

3 Cone— shaped cluster CALOl 

In the case of calorimeter event selection the detector does not discriminate 
the particles moving near the same direction, i.e. producing a cluster. In the 
first order of perturbation theory the cluster can be formed only by two particles: 
photon and electron (or photon and positron). For the definiteness we shall talk 
about the positron cluster, although the same will be refer to the electron one. 

If a photon and positron belong to a cluster (7, e + ) G CL, the detector mea- 
sures the total energy of the cluster, and its position is determined by the 
positron position. If a photon and positron do not form a cluster (7, e + ) CL, 
the event looks exactly the same as in the case of BARE1, i.e. there exists 
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the detection threshold on the positron energy and the photon is entirely not 
detected. The cutoff on the positron energy can be written symbolically in the 
form 

1 X c 1 X c 

J dx + J [(7, e+) e CL]dx ee J dx- J [(7, e+) g CL]dx . (14) 

x c 

It is convenient to use the left side of Eq.(14) in order to take into account 
the initial-state emission and right side of Eq.(14) for the description of the 
final-state one. As follows from Eq.(14), in the calorimeter event selection the 
correction to the cross-section can be represented in the form of two-term sum. 
One of these terms is a universal, i.e., it does not depend on the cluster specific 
and another one is determined by its form and size 

Si = s7 + s c ; . (15) 

The contribution of the universal part S^™ can be obtained with the help 
of the formula (13). For this one must retain without change the correction 
caused by the initial-state emission and integrate the correction, caused by the 
final-state emission, for x going from to 1. The result has the following form 

P% 1 

sr = f / ^{-Ag+ / [((i-i)p 1 ( :C )+i-x+i^-L 2 )A<?+(5r ) -or ) )£ 3 ]rfx+ 

Z7T J Z \ 1 — X 

+ / [(!-* + V^ Z i) A 42+ 1 -7^(B A -9 2 )L A ]dx\+ 
J 1 — x 1 — X > 

Pa 1 

+V- I + I [« L " + 1 - x + \ t ^-L 2 )8i x) + 9 ( * ) L 3 ]dx+ 

ZTT j „ Z ^ J ' 1 — X 

p\ Xc 

+ / (1-* + -^— Li)dz}, (16) 
I x 

where Li are obtained from Li by the substitution p 3 instead of and 1 instead 
of p2- In the case of symmetrical ring-shaped detectors it is necessary to put 
Pi = Pi — 1 on the right side of Eq.(16). 
We see that owing to the well-known relation 

1 
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the final-state emission does not give rise to large logarithm L in the cross- 
section in accordance with the Lee-Nauenberg theorem about the cancellation 
of the mass singularities when summing over all possible states [21]. 

In order to obtain the part of the correction, determined by the shape and 
size of the cluster, one may use the simplified form of the bremsstrahlung differ- 
ential cross-section, that describes the so-called semicollinear kinematics (see, 
for example, Refs. [14,22]). In the case of single-photon emission it simply 
means, that one can neglect the electron mass everywhere. The corresponding 
expression for the dimensionless quantity S, related to the photon emission by 
positron, is determined as follows 

^ adipdzdzi(l + x 2 ) . 1 X ]d 

4kQ 2 z(z\ — XZ) Z\ + Z + 2yf Z\ZCOS(f Z\ + X 2 Z + 2XyfZ\ZCOSLp 

(17) 

The cutoff on the z, Z\ and if variables depends, of course, on the cluster 
shape. The CALOl cluster is the cone with the cone semi-angle 6 = 0.01 around 
the final-positron 3-momentum. In order to include the initial-state emission it 
is necessary to have (according to Eq. (14)) the cutoff on these variables when 
photon and positron belong to a cluster 

i j— I- \ n A 2 x 2 (l - x) 2 - (Jz[- xJz) 2 x 6 
\Jz\-xyJz < x{\-x)\ -1 < cosif < -H '- — ^—t-i A = — 

lXsJZZ\ V\ 

(18) 

For the calculation of the correction caused by the final-state emission it is 
necessary to use the cutoff when photon and positron do not form a cluster 

— 1 < cos ip < 1, \\[z[ — x\fz\ > Xx(l — x); 
„ X 2 x 2 (l — x) 2 — [ x fz{ — xJ~z) 2 . , 

1>C0S<^>-1 + '- y — L , \y/zi-Xy/z\<\x(l-x). (19 

ZX-yJZZl 

We pay attention to the fact, that formally the cuts (18) and (19) determine 
also the collinear and semicollinear kinematic regions, respectively (in the case 
of the final-state emission [12,13,22]). The only, but very essential difference is 
that, in the considered case , the parameter A is of the order of unity, whereas 
when determining the kinematic regions, it is required that this parameter was 
much less than unity. 
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Integration of Eq.(17), as well as a similar formula including photon emission 
by an electron, over regions (18) and (19) leads to the following contribution to 
the first-order correction 



Oi r 1 -|- X A r dZ r 

= — / dx / d Zl (^ + ^)^(z h z;X,x), 

2tt J 1 — x J z l J 



dz 



i/ ' 



(20) 



a r 1+ X 
lf ~ 2TT \ 1-x 



dx 



{ / So* 



xp\ — z 



P3 



+ M + j ^(ln 

a 2 Z 



X — Z 



1 - Z 



+ 1- + 



+ 



/ ^[0(ag-,)-0(,-62)]Z 4 + / "(in 

1 ^2 Z 



xpi — z 



pl-z 



P2 



+i- + 



J z z J > z =p zxX^z — X z x(l — x) 



In 



(21) 



where the integration limits over z (in the squared brackets) and Z\ (in the round 
brackets) are defined with the help of \I> and ^ quantities 

\fr = [a 2 ,al}(x 2 z + ,x 2 ) + [b 2 ,a 2 }(x 2 z + ,x 2 Z-) + [bl,b 2 }(x 2 pl,x 2 Z-), z± = (y/z±X(l-x)) 2 , 

^ = [d 2 ,dl](x 2 z + ,x 2 pl) + [6 2 , a 2 ](x 2 z + , x 2 Z-) + 6 2 ](x 2 /0 2 , x 2 z_) , 

z± = (>fz ± A(l - a:)) 2 , ao = p2 , 6o = ^4, a = max(p2, 1 + A(l - x)) , 

6 = min(p^ , /?3 — A(l — x)) , ao = max(l, /?2 — A(l — x)) , a = P2 + A(l — x) , 

6 = /?4 — A(l — x) , 60 = min{p^ + A(l — x) , p%). (22) 

The $ and F functions beloning to the right side of Eq.(21) are determined as 
follows: 

z — Z\ 



$ = 



F= 2 - 



1 



+ 



1 



7T \Z\ — XZ Z — Z\ 
1 1 



arctan <^ 



IT \Z\ — XZ Zi — X"Z 

R = 



arctan 



(JZ[- Xjz) 2 



R 
R 



Zi — X A Z 



A 2 x 2 (l — x) 2 — (y/zi — x-Jz) 2 



\ (v^r + x^~zf - a 2 .x 2 (i - .x) 2 ' (23) 

The formula (21) is valid in the general case of the asymmetrical detectors . 
In the symmetrical case of the wide detectors it is necessary to put p± = p% and 
p 2 = 1 in Eq.(21) 
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4 Pyramidal cluster CAL02 



The cluster CAL02 has a pyramidal form. Its section by a plane perpendicu- 
lar to the beam axis represents a rectangle in the angular coordinates (0, ip), and 
the centre of it is determined by the position of the positron 3-momentum. The 
length of this rectangle along polar axis is equal to 20q, and along an azimutal 
one - 2$. Their values are 

O = 0.051/16, $ = 3tt/32. 

The cuts, for the cases when a photon and positron form or do not form 
cluster CAL02, are very simply formulated in the terms of the positron and 
photon variables, namely 



\0+ - 7 | < 0o j 1 > cos ^7 > cos $ 5 (7 j e+ ) £ CL, 



(24) 



and also 

\0 + — 7 | > 0o, 1 > cos<£> 7 > —1; 
1 0o - 7 | < 0o, -1 < cos<^ 7 < cos$, (7,e + )^CL, (25) 

— * . 

where 7 = k / (e(l — x)), and (p 1 is the angle between two-dimensional vectors 
k 1 - and 

It is necessary to express the conditions (24) and (25) in terms of electron and 
positron variables z, z\ and <p. Note in this case, that the law of the transverse 
momentum conservation leads to the following relations: 

0_ cos ip + X0 + 



cos Lp 1 = 



7 = —^—J0 2 _ + x 2 0% + 2x0+0. cos (f . (26) 
1 — X 



(1 - x)0 1 ' 

Using these relations we take for the case when photon and positron form a 
cluster 



A < if < B, A = max 



-1 = 



(1 - xf{0 + - 0o? - (0- - x0 + Y 



B = min 



(l-x) 2 (0 + + e o ) 2 -(0^-x0 + y 
2x0+0- 



2x0+0- 



cos($-£) 



-1 



sin S = 



X0 



+ 



\ 0- 



sin 
(27) 
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From Eq.(27) it follows that there are three different kinematic regions if 
photon and positron belong to the cluster CAL02. In the first region we have 

(1 - xfiJzl - xX) 2 - x 2 Jz - Jz{) 2 

-1 + v — - — y y - J± - <cos(p<- cos($ - 5) , 

2x l yfzz\ 

x 2 z+ < zi < x 2 J\. (28) 

In the second one 

- 1 < cos if < - cos(<i> - 5) , x 2 z + > zi> x 2 J 2 , (29) 

and in the third region 

(1-s) 2 (^+iA) 2 -i 2 ^-^) 2 

— H - > cos ip > - 1 , x z + < z\ < x J , , 

2x l y/zzi 

(30) 

where 



J (±) = ^ 2 (1 - x) 2 X 2 sin 2 $ ± (l-x)A(l-2a:sin 2 2)] , 

= 1 - 4ar(l - ar) sin 2 2 , A = ^ . 

If the photon and positron do not form the cluster CAL02, then it is evident 
that 

- 1 < cosip < 1 , zi> x 2 j\ , zi < x 2 z_, (31) 

and in the opposite case 
^ (1 - x) 2 (^/zj- xX) 2 - x 2 ^ - ^/zj) 2 

-1 + - ' Y — >COStf>-l, -COs($-5) < COS(Z? < 1, 

x 2 z + < z\ < x 2 j\ , (32) 

- cos($ - £) < cos<^ < 1 , £ 2 z + > z\ > x 2 J 2 , (33) 

, , (1 - x) 2 (yfzT + zA) 2 - a: 2 ^ - ^z[) 2 2 2 2 

- 1 H < COS <y2 < 1 , X Z < Z\ < X J . 

2x A yfzzi 

(34) 

Integration of Eq.(17) leads to the following result for the contribution to the 
part of the first order correction which depends on the shape and size of the 
cluster CAL02 



/ \^-dx \% /dzi[(^i + *i)$i + (^ 2 + ^ 2 )$ 2 + 
J Q 1 — x J z z J 
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y — 

^ ~ 27T 



+ ($3 + $ 3 )$ 3 ]- 



7T \Zi — XZ Z — Z\ 
Si; = ^ / Y^dx{j ^ / dz 1 [(^ 1 + f OFx + (tt 4 + ^2 + 



(35) 



~ 2/1 
+(^ 3 + ^ 3 )F 3 ]- 



1 



7T \Z\ — XZ Z\ — X^Z 



+ 



Pi 



dz 



+ / ^[#(maa;(l, Z2 ^ — z) — 9(z — min((/?4 + A(l — x)) 1 , Pz)]L,4+ 



2 2n 



+ 



1 dz 



XP2 — Z 



P2 - Z 



v(-) 



+ i- + 



/ 



max(l,Z2 ^) 



Z^ 



In 



(^P4 - z)(J+ - z) 



(zJ 2 - z){p\ - Z) 



+ 



min(pl,z { 3 } ) ^ 



Pi 



f 2 
^ 3 



(^J 2 - Z)(p| - z) 



/ 



Z — X 



1 - z 



+ i_ , (36) 



where 



4 ±} = (pi ± (1 - ^)^) 2 - 4a: (1 - a;)pi(pi ± A) sin 2 5. 



The $i and functions, entering in Eqs.(35) and (36), are defined as follows 

>(-) 



3>i = arctan Q- — arctan 77, <3>2 = arctan^ \ $3 = arctan 



Q 



(+)' 



iq = arctan ) ^ , F2 = arctan C, F3 = arctan ) , , 

0/ 0/' 

Hi 



rj = Ri cot 



C = Rf cot 



-(±) = ^/q(±) r . = (v^Z v 7 ^) 2 _ (y^I-^v^) 2 

2 



z — Zi 



/ 



zi — arz 



(±) _ 



^7 



a; 2 (v/z + V ^I) 2 - (1 - x) 2 (yfz{ ± ^) 2 



(37) 



^ (1 - xfiy/zi ± x\y -x\j~z- yii) 2 • 

The and ^ quantities are inserted in Eqs.(35) and (36) in order to indicate 
the integration limits over z and z\ 

*i = [min(plzt ] ) , p 2 2 ](x 2 J 2 + ,x 2 z + ) + [6 2 , z^ } ](^, x 2 z+) , 

^2 = [max(p\, z[ + ^) , /^(a^z^ £ 2 ) + [fr 2 , max (pi, z[ + ^](x 2 z + , x 2 J 2 ) + [ / o 2 , fr 2 ]^ 2 /? 2 ;, a; 2 J 2 ) , 
^ 3 = [a 2 , max(p\, z[ + ^)](x 2 J 2 , x 2 ) + [p 2 , a 2 ](a? 2 J 2 , £ 2 z_) , 
^4 = [max(pl, z{ + ^), pl](x 2 J+,x 2 ) + [min(/)], zj"'), raas^, ^ + '](i 2 4i 

+ [/9imm( / 9iz^ ) )](a: 2 ^,x 2 J 2 ) , 
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^i = [(pi - X(l - x)) 2 ,max(l,z { 2 } ](x 2 j\, x 2 pl) + [z\ \ a ](x 2 j\, x 2 z + ) + 

[(p 4 -X(l-x)) 2 ,z ( f ) )](x 2 plx 2 z + ) , 

^ 2 = [z { 2 + \a 2 ](x 2 z + ,x 2 p 2 ) + [(p 4 - A(l - x)) 2 ,zi +) )](x 2 z + ,x 2 J 2 _)+ 

+ [min(pi ,(p 4 -A(l-z)) 2 ](z 2 p 2 ,z 2 J 2 ) > 

^3 = + A(l - x)) 2 ,4 +) ](^V-,^ 2 P2) + [6o, (P 2 + A(l - x)) 2 ](x 2 J 2 + ,x 2 z + ) , 

^ 4 = [zi + \max(l,zi~ ) )](x 2 Jl,x 2 pl) + [4 _) , 4 +) ](^4^ 2 J -)+ 

+ [min(p3, 4 +) ) 5 4 -) ](£ 2 /C>4,£ 2 J 2 ) • 

It is necessary to replace A by A in the used here quantities a, 6, ao, &o- As before, 
in the case of symmetrical wide detectors, in Eqs.(35) and (36) one must replace 
P4 by pz and pi by 1 and simultaneously put to zero those terms where the upper 
integration limit over z becomes less than bottom one. 

5 Leading second— and third— order corrections 

As we have seen in the previous Sections, the large logarithm L enters only in 
the universal part of the radiative correction. Such situation takes place also for 
leading contributions in the higher orders of the perturbation theory. Therefore, 
the analytical formulae for high-order leading radiative corrections are universal. 
They are applicable for both CALOl and CAL02 event selections. 

In order to obtain the leading contributions in the second and third orders 
we use the method based on the electron structure functions in the singlet and 
nonsinglet channels [16] and effective coupling constant [7] 

D{x, oteff) = D NS {x, oteff) + D s {x, a eff ) , (38) 

where the effective coupling constant is defined as the integral of the running 
electromagnetic constant in one-loop approximation 



15 



The iterative form of the nonsinglet component of the structure function can 
be presented as 



D NS {x,a eff )=5{l-x)+Y.ii 



1 fOL e ff 



i k\ V 2tt 



Pi(x)<8» •••<8»Pi(x) = Pi(a;)®* , P 1 (x)®P 1 (x)= f P 1 (t)P 1 (^)^ . (40) 
v . ' % u/ r 

The singlet component of the electron structure function including third- 
order contribution can be written as follows 



(41) 



R{x) = 



1 — x 
3x 



(4 + 7a: + 4a; 2 ) + 2(1 + x) In a; . 



The nonsinglet component of the electron structure function describes the 
photon emission and electron-positron pair production neglecting the final- 
electron identity, whereas the singlet one just responds for the identity effect. 

With the required accuracy the electron structure function reads 



aL 



1 /<xL\2 



D(x, L) = 6(1 -x) + —P^x) + - {^) G(x) 



1 /aL\3 



3V2tt 



G(x) = P 2 (x) + ^P 1 (x)+R(x), F(x) = ^P 3 (x)+P 2 (x) + ^P 1 (x) + ^R(x)+R x (x), 

Pi(x)=P®\ R{x) = P l ®R{x). (42) 

The functions Pi(x), P%(x) and P P (ir), entering in the right side of Eq.(42,) can 
be written as follows [7] 

R P (x) = ll + 2 m (! " x)]R(x) + (1 + x) [4L i2 (l -x)- In 2 x] + i(-9 - 3a:+ 

2 3 

+8x 2 ) lnx + -( 8 + 8a: + 3a? 2 ) , 

3 x 

Pi(x) = Ai6(l - x) + Qi(x)e(l-x -A), z = 2, 3, 
A 2 = (2 In A + ^) 2 - 4C 2 , A 3 = (2 In A + ^) 3 - 12( 2 (2 In A + \) + 16Cs , 



e 2 (x) = 2 



6 3 (a;) = 12 



1 + a; 2 
1 — x 



In 



1 + x z 
1 — X 
1 — X 



(In- — + -) + -(1 + re) In a: - 1 + x 

0C Zi Zi 



ln(l-a;) + -ln 2 a:+^ln^ — +^"C 2 -3(l-x)[l+ 

x 6 4 a; 16 
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3 rl (1 - x) 4 
+41n(l - x)] + -(5 - 3rr)ln£ + 6(1 + x) -In a; In- — + L i2 (l - x) 

The functions Pj(x) satisfy the condition 

l 

J Pi(x)dx = , 
o 

which is nothing but the statement of the Lee-Nauenberg theorem [21] for the 
leading radiative corrections in the singlet channel in terms of the electron struc- 
ture functions. 

The factorization form of the differential cross-section in the frame work 
of the impact representation [23] determines the leading contribution to the 
radiative correction in the case of calorimeter event selection in the following 
form 

OO 7 1 1 

E L = J ^ J d Xl J dx 2 D(x u L)D(x 2 ,L) A^A^ . (43) 

x c x c jx\ 

The integrand in Eq.(43) includes only those ^-functions which depend on x\ 
and x 2 . It corresponds to account for the initial-state emission. The final-state 
emission does not produce the leading contribution according to the mentioned 
above Lee-Nauenberg theorem. 

Combining Eqs.(42) and (43) we obtain the leading contribution in the second 
and third orders 

2 oo 7 1 1 

£2 = £2 / f L *{ I [(Af 2 ) A 31 +Ag ) A 42 )G(x)+2 / dx 1 P 1 (x)P 1 (x 1 )A^A^ 

x c Xc jx 

(44) 

3 00 7 l-i 

^ = h I f L 1 / [^A^Aa! + AWA 42 )FW + 

x c 

+ i / dx 1 P 1 (x)G(x 1 )(A$A$ ) + A^A^]dx) . (45) 

x c /x 

In order to eliminate the 6 - functions, entering in (44) and (45), and place 
specific limits, the additional work is needed. Omitting all intermediate calcu- 
lations let us write the second-order leading correction as follows 

S 2 = S 2R + S 2(+) + S 7 , (46) 
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where the first term on the rihgt side of Eq.(46) is caused due to the emission 
of the real and virtual photons as well as electron-positron pair production by 
an electron, the second term - by the positron, and the third one - due to the 
emission of one photon by an electron and one photon by a positron: 



n 

l a 2 r r dz 
^ = 8^ 



dz 



p\ 



m 23 



P3 



(47) 



Pi 



pI 



y 7 _ _ 
7 4tt 2 



mi4 



{- / S^w- / S --)]- / 

„2 ^ rni.1 ra rr, 



PA 



Pi 



+ / -^B 



m 12 



P2 



, (48) 



dZ^cy^ 

—Z-LF n 



Pi 



P2 



d^ T 2 F (V± 

2 9\ „ ' ^ c 



mi2 



P2 



pI 



x c ) + 



1/ 



X c p2 



+ C[y/z 



Pa 



pI 



XcP3P2 



P2 



dz r0 

z z 



y/z x c p 3 



Pa 



+ c 



Pa 



Z \ z 



XcPA 

Z X c p 2 



F n 



P'd Pa 
\fz x, 



+ 



9\ „ ' 



Pa ' 



+ 



P3 Vz ' ^ P2 Pi 



Z \ z 



(49) 



where 



rriij = max(xlp 2 j, p 2 ) , A = — J (P2(x)+-Pi(x))dx+ J R(x)dx, B(y) = J G(x)dx. 



F 2 (x) = J P 2 (y)dy = -2x-^-+(x+^-) In 



ar 



X 



4 -V-. 2 ^" (1 _,)4+ 41n ( 1 -) ln T^+ 4 ^) 

F,(x,y) = F,(x)-F,(y), 



F g (x) = | Pi(x)p cfe = -^+(2x+x 2 )lnx+^ c + yj In 



(1-xf 



2x r +°^ 



-2x - y) ln(l - x c ) + 4L z2 (» + 4L l2 ^ — , x c < x < 1 



3/ /■ 



g(y) - g 



tjj ( 



X 



/X 
P\{x)dx = x + — + 21n(l — x\ x < 1 



(50) 
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By analogy with Eq.(46) the third-order leading correction can be written as 
follows 

£3 = s;+s;+e 2 1 + s; , (51) 

where the upper index shows the number of real and virtual particles emitted 
by an electron and the bottom one - by a positron. For the case when three 
additional particles are emitted by one of the fermions, the corresponding con- 
tribution to £3 reads 



(\ 3 Pi ,7 x c 1 

tH { / ^ 3 [" 2 / F p (x)dx + 2 J F r (x)dx 

' pi Z Xc 





p\ ^ Vz/P3 Pi ^ ^ V^/Pi p\ ^ z Vz/p 2 



— J —^L 3 J F c (x)dx— J —^L 3 J F c (x)dx+ J -^L 3 J F c {x)dx^ 

X c 

(52) 



mi4 



m.12 



F P (x) = \p,{x) + l -P 2 {x) + A Pl(x ) } 

F r (x) = 1r(x) + l -R\x) , F c (x) = l -F(x) . 
In the case when both fermions radiate simultaneously, we have 



£1 + E; = (|-) 3 { / ^L 3 [- J (p 3 (x) + ?P 2 (aO) dx + / 

p2 ' x c 



pl dz 



y — I/ 3 J H(x,x c )dx— J —L> 3 J H(x,x c )dx+ J —1^ J H(x, 



m 23 



mi4 



m 12 



+ 



y ^-^ 3 y dxP\{x)N I J + (p 3 <-> p 4 ) 



X C P3P4 



XcPi/^/z 



\X /04 



+ 



+ y ^r 3 [ y 7 dxP^N 

X ^ Z x c p 2 /^ 



X c \l z 



X p 2 



dx + (p 2 <-> 1] 



- J L 3 [ J dxP 1 (x)Nl^y^\dx + (p 3 ~p 2 ) 

X ^P 3 x c p 2 /^~z \ X P2 ' 

- J %L 3 \ "7 dxP^N (^y^-)dx + (p 4 ^l)]}, (53) 



X c p 4 /yfz 
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where 
H(x,x c ) = Pi(x) 



+ 3 5 



x J 2 



1 

+ 2 S 



tjj f 



X 



(P 2 (x) + R(x)), 



f(y) 



N(x; y) = N(x) - N(y), N(y) = + ±g(y) + ir(j/), 

-F 2 (y), r(y) = J R(x)dx = -— + y + ?/ 2 + g?/ 3 - ^ + 2y + In 



When writing the formulae of this Section, we represent the restrictions on the 
angles and energies of the detected particles with the help of definite integrals 
using such relations as, for example, 



j 9^6^ 6^ dz dx dx\ = J dz J dx J dx\ 

XcP3 x c p 3 /^/z x c /x 



(54) 



It is necessary to bear in mind that the upper limit of the integration over the 
variable z always must be greater than bottom one. In opposite case the integral 
must be putted to zero. 



6 Numerical results 

At LEP1 conditions the limiting angles of the ring-shaped detectors have 
different values for various versions of event selection. In the cases of BARE1 
and CALOl the wide detector has 

0i = 0.024 , 6> 3 = 0.058 , 

and narrow one - 

9 2 = 0! + h , 4 = 03 - h , h = . 

o 

The wide detector (ww) for the event selection CAL02 coinsides with the nar- 
row one (nn) for the case BARE1, whereas the limiting angles of the narrow 
detector are determined as follows 



6> 2 = 9i + 2h , 4 = 3 - 4/i . 

The Born cross-section, determined by the formula (4), is equal to 
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o B = 175.922 nb for ww BARE1 CAL01, 

o B = 139.971 nb for ww CAL02, nn BARE1 and CAL02, 

o B = 103.299 nb nn CAL02. 

The results of our calculations of the radiative QED corrections are presented 
in the Tables I— III, where the vacuum polarization is switched off. For the 
comparison we also present numbers in the Tables I and III, that are obtained 
with the help of the MC-generator BHLUMI using the exponentiated Yennie- 
Frautchi-Suura factor for the higher-order corrections. 



x c 


bhlumi ww 


WW 


nn 


wn 


calol 


0.1 


166.329 


166.285 


131.032 


134.270 


0.3 


166.049 


166.006 


130.833 


134.036 


0.5 


165.287 


165.244 


130.416 


133.466 


0.7 


161.794 


161.749 


128.044 


130.542 


0.9 


149.925 


149.866 


118.822 


120.038 


calo2 


0.1 


131.032 


130.997 


94.666 


98.354 


0.3 


130.739 


130.705 


94.491 


98.127 


0.5 


130.176 


130.141 


94.177 


97.720 


0.7 


127.528 


127.491 


92.981 


95.874 


0.9 


117.541 


117.491 


86.303 


87.696 



Table I.SABH cross-section at LEP1 conditions with 
first- order correction. Vacuum polarization 
is switched off 



From the Table I one can see, that at LEP1 conditions, the SABH cross- 
section including the first-order correction, obtained by the means of BHLUMI 
generator systematically exceeds our results in all interval of values for parame- 
ter x c . The relative difference amounts to about 0.03%. The possible reason of 
this defference, as it was noted in Section 2, is that we neglect the contributions 
proportional to 9\ when calculating the radiative corrections. Within the accu- 
racy of 0.1% the contribution of the omitted by us terms is insignificant. But it 
requires the additional systematic investigations if the accuracy will be 0.05% 
(and it had actually been achieved in the laboratory L3 [24]). 
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The absolute values of the second- and third-order leading corrections (in 
nb are presented in the Table II. The second-order correction is divided into 
the contributions related to the electron-positron pair production and double- 
photon emission. One can see that photonic corrections dominate in the second 
order of the perturbation theory. The third-order correction includes the con- 
tributions due to three-photon emission and pair production accompanied by 
single-photon emission. 





calol 




calo2 




x c 


WW 


nn 


wn 


WW 


nn 


wn 


correction due to pair productiion 


0.1 


-0.046 


-0.045 


-0.024 


-0.045 


-0.047 


-0.024 


0.3 


-0.046 


-0.045 


-0.024 


-0.045 


-0.047 


-0.024 


0.5 


-0.048 


-0.046 


-0.025 


-0.046 


-0.047 


-0.024 


0.7 


-0.069 


-0.059 


-0.042 


-0.059 


-0.051 


-0.036 


0.9 


-0.137 


-0.111 


-0.102 


-0.111 


-0.085 


-0.075 


second order photonic correction 


0.1 


0.788 


0.708 


0.302 


0.708 


0.668 


0.255 


0.3 


0.680 


0.634 


0.195 


0.634 


0.627 


0.187 


0.5 


0.474 


0.487 


0.005 


0.487 


0.554 


0.073 


0.7 


0.293 


0.317 


-0.164 


0.317 


0.396 


-0.092 


0.9 


0.866 


0.738 


0.373 


0.738 


0.628 


0.271 


third order correction 


0.1 


-0.041 


-0.036 


-0.002 


-0.036 


-0.034 


-0.001 


0.3 


-0.046 


-0.040 


-0.007 


-0.040 


-0.037 


-0.003 


0.5 


-0.044 


-0.039 


-0.006 


-0.039 


-0.037 


-0.005 


0.7 


-0.023 


-0.022 


0.012 


-0.022 


-0.027 


0.008 


0.9 


0.021 


0.013 


0.049 


0.013 


0.002 


0.038 



Table II Absolute values of second- and third-order 



corrections to SABH cross-section at LEP1 conditions (in nb) 



In the Table III we present the total SABH cross-section at LEP1 with ac- 
count for all corrections calculated in this paper and also the results of the 
corresponding calculations made with the help of the MC-generator BHLUMI. 

As to comparison of our calculations and MC-generator BHLUMI ones for 
the radiative corrections in the second and third orders, it is necessary to note 
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that BHLUMI results are based on the exponentiated form of electron structure 
function whereas our ones - on the iterative form. The corresponding effect due 
to these different forms increases at large values of parameter x c as one can see 
from the Table III. 

For more efficient comparison, it is necessary to have either the analytical 
calculations with the exponential form of electron structure functions or MC- 
calculations that do not use the exponentiation. For MC generator BHLUMI 
such calculations exist in the case of BARE1 event selection [3]. The correspond- 
ing comparison have been recently done and the agreement has been obtained 
at the very high level [25]. 





calol 


calo2 




x c 


ww nn 


wn ww nn 


wn 


total cross— section (analytical calculation) 


0.1 


166.968 131.659 


134.546 131.624 95.253 


98.584 


0.3 


166.594 131.382 


134.200 131.254 95.032 


98.285 


0.5 


165.626 130.818 


133.440 130.543 94.647 


97.764 


0.7 


161.950 128.280 


130.348 127.727 93.299 


95.758 


0.9 


150.616 119.462 


120.358 118.131 86.848 


87.930 


total cross— section (bhlumi— generator) 


0.1 


167.203 


131.835 95.458 


98.834 


0.3 


166.795 


131.450 95.233 


98.539 


0.5 


165.830 


130.727 94.841 


98.020 


0.7 


162.237 


127.969 93.520 


96.054 


0.9 


151.270 


118.792 87.359 


88.554 



Table III Total SABH cross-section at LEP1 conditions 
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